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Abstract 



We study WeyssenhofF spin fluids in Bianchi type-I cosmological models, 
, within the framework of torsional /(i?)-gravity; the resulting field equations 

■ are derived and discussed in both Jordan and Einstein frames, clarifying the 

role played by the spin and the non-linearity of the gravitational Lagrangian 

^ f{R) in generating the torsional dynamical contributions. The general conser- 

vation laws holding for /(i?)-gravity with torsion are employed to provide the 
^ ■ conditions needed to ensure the preservation of the Hamiltonian constraint 

t^l^ and the consequent correct formulation of the associated initial value problem. 

Examples are eventually given. 

^ ■ Keywords: Spin fluids, /(i?)-cosmology, torsion 
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^ ■ 1 Introduction 

■ Despite most of the success General Relativity (GR) enjoyed, there are some prob- 
lems, at both ultraviolet and infrared scales, that GR is unlikely to solve, and 
therefore any of its enlargement is welcome: because the structure of the Einstein 
gravity is based on the assumption of least-order derivative field equations, by re- 
laxing it means going to higher-order derivative models; these extensions however 

I suffer from some arbitrariness, as nothing tells us what is the order derivative of 

the field equations. 

However, among the possible higher-order derivative field equations, or actions, 
there are two that are somehow special: one of them is based on the assumption 
of conformal invariance, and its beauty comes from the fact that this simple prin- 
ciple of symmetry is enough to severely restrict the possible actions to just a few; 
the other is based on the assumption of considering all gravitational information 
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contained in the Ricci scalar curvature alone, then taking the most general form 
f{R) for the gravitational action. Although the principle of conformal invariance 
is a beautiful and simple symmetry principle, nevertheless the universe does not 
display such a symmetry and a mechanism of spontaneous symmetry breaking 
must be invoked, and without it, it is preferable to investigate f{R) gravity as a 
candidate to solve the problems we have mentioned above. 

In fact, f{R) gravitation has been so far useful to address astrophysical and cos- 
mological problems that we could otherwise faced only with the introduction of 
exotic form of matter and energy, the Dark Matter and Dark Energy; these forms 
of matter however are still a theoretical mystery. By converse, f{R) gravitation 
has been able to account for both cosmic speed-up and missing matter without 
the introduction of new elements, as it is discussed in [ll[2l[3]. 
On the other hand, the structure of the geometrical background must also take 
into account the matter fields we find in it: matter fields, according to the Wigner 
theoretical classification in terms of mass and spin, or the heuristic fact that they 
possess translational as well as rotational degrees of freedom, must have both en- 
ergy and spin density; in a framework in which, as it appears to be reasonable, all 
physical quantities are coupled to geometrical ones, the energy is to be coupled 
to the curvature in the same spirit in which the spin is to be coupled to torsion. 
Thus we retain torsion not be neglected, likewise we would not neglect curvature, 
coupling the former to the spin, and the latter to the energy, and giving rise respec- 
tively to the Sciama-Kibble field equations, beside the Einstein field equations, as 
it has been discussed in [4J. 

Now according to this imprinting, a theory of the f{R) type in which the Ricci 
scalar curvature R contains both metric and torsional degrees of freedom is the 
most straightforward extension we can take into account; an important conse- 
quence of the non-linearity of the function f{R) is that we have non- vanishing 
torsion even without spin in general (that is if the energy-momentum trace is not 
constant O El [71 [8]). Then torsion may give rise to a repulsion ensuring singularity- 
free and accelerated cosmological models ([9]). 

Besides, in a recent paper |10) we have studied in detail this model in presence 
of Dirac fields, in order to better exploit the coupling between torsion and spin 
within the framework of /(i?)-theories finding that, as it might have been some- 
what expected, the angular degrees of freedom carried by the spin are amplified 
by its coupling to torsion in such a way that one is forced to deal at least with 
axial configurations in order to avoid consistency problems. Thus in another work 
[llj we have focused on the Dirac field exploring cosmological applications for the 
simplest anisotropic universe given by the Bianchi type-I (BI). 
In the present work, we explore yet another situation, in which what is supposed to 
constitute the macroscopic approximation of the Dirac quantum field, the Weyssen- 
hoff spin fluid, is studied in the context of BI universes. 
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2 /(i?)-gravity with torsion: generalities 



The dynamical gravitational variables of /(i?)-gravity with torsion consist in pairs 
((7, r), where g is a metric tensor and F is a metric-compatible linear connection. 
The field equations of the theory turn out to be O El [3 [5] 



f{R)R^J - -f{R)gij = S^, (2.1a) 



rp h 



(2.1b) 



where Rij and T^^^ are the Ricci and Torsion tensors associated with g and T, 
and Sjj and S'j are the energy- momentum and spin density tensors of the matter 
fields. From (f2Tbl) . it is seen that we can have non-vanishing torsion even in 
absence of spin density. Following previous works [51 El El El HD] , we suppose that 
the trace of equations ()2.1ap 

f'{R)R-2f{R) = T. (2.2) 

gives rise to an invertible relation between the Ricci scalar curvature R and the 
trace S of the energy- momentum tensor. Also, we assume that f{R) 7^ kR? 
(the case f{R) = kR? is only compatible with the condition S = 0). Under the 
assumed conditions, from equation (|2.2p it is possible to derive the expression of 
R suitable function of S, namely 

R = F{T.) (2.3) 

Making use of the Bianchi identities, it is possible to derive the general conservation 
laws 

ViS^^' + T,T}^ - ^piT^P' - ^SsuR'''' = (2.4a) 

VhS'^^ + ThS'^^ + T:'^ - S^' = (2.4b) 

under which the energy-momentum and spin density tensors of the matter fields 
must undergo once the matter field equations are assigned (TOj- In equations 
(j2.4p the symbols Vj denote covariant derivative with respect to the dynamical 
connection F while i?***-' is the curvature tensor of F. Indices are raised and 
lowered by the metric g. 

The dynamical connection can be expressed as 

F,/ = f>-K,> (2.5) 

where T-^ denotes the Levi-Civita connection induced by the metric gij and K-^ 
indicates the contortion tensor 1121 



^i/ = ^(-^^/ + ^/i-^%) (2-6) 
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Using the second set of field equations ()2.4bp , we obtain tlien the following repre- 
sentations 

i^,/ = K,> + 5,> (2.7a) 

^"■^^ 2/'(i?(S)) {~^'^^ ~ ^^'■^) (2.7b) 

Kij'':=-f,6'^ + fpgP^gij (2.7c) 



which, together with equation ()2.5p . turn out to be very useful in order to separate 
the purely metric terms from the torsional ones in the field equations ()2.ip . 



3 Coupling to spin fluids: fleld equations and conser- 
vation laws 

We consider a Weyssenhoff spin fluid, characterized by an energy-momentum ten- 
sor of the form 

J^ij ^ Jjipj + p (^Uijjj _ gij-J la) 

and a spin density tensor given by 

= S.jU'' (3.1b) 

where W {U^Ui = 1) and denote respectively the 4-velocity and the 4-vector 
density of energy-momentum, while Sij the spin density of the fluid (see, for ex- 
ample, [13l[ll] and references therein). The 4-velocity and the spin density satisfy 
the convective condition 

SijW = (3.2) 

It is easily seen that the relations (|3.2p . together with equation (|2.7bp . imply the 
identities 

5/'^ = -5'>^ = (3.3) 

Making use of ([331) as well as of equations ([23]), ([Ml), ([221), (l3lbj) and ([32]), we 
can express the contracted curvature and the scalar curvature respectively as 

Rij = Rij - 2Vjfi - Vhf^Qij + 2fifj - 2fhf^gij - j,fhS^jUi 

N 1 (3.4a) 

and 

R = R- 6V,r - 6fif + -^Ij^SP'^Spg (3.4b) 
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where now Tj = -fjr^- In view of this, substituting equations (|3.4|) in equations 
()2.ip . we obtain Einstein-hke equations of the form 

~ 1 P _ 1 ^ 1 f 3 d(p d(p fi dip d(p f^f^ 

-^V^'^Qij - V{^)gi^ + ^T)^5^C/, + {-SjiU'' + % - S'^jU^ 

1 1 

(3.5) 

where, using equation (|2.3p . the scalar function 

^ ■■= nn^)) (3.6) 

and the effective potential [S] 

V{^) := i [^F'\{frH^)) + ^HfTH^)] (3.7) 

have been introduced. It is a straightforward matter to verify that the conservation 
laws for the spin ()2.4b[) amount to the validity of the anti-symmetrized part of the 
Einstein-like equations ()3.5p 



V;^ + nS''^U^ - fhS\Uj + S,,- - = (3.8) 



Saturating equations (j3.8p with [/* we obtain the explicit expression for the 4- 
vector density of energy-momentum 



P, =pU,- ThS'^j - Vh [Sij U^) W (3.9) 

where p := U^Pi. inserting equations (|3.9p into (|3.1ap we get the explicit form of 
the energy-momentum tensor 



^ij = {P + P) UiU, -pgij - UiThS^j - UiVh [SkjU'^j (3.10) 
The significant part of the Einstein-like equations is then the symmetric one 

R _if?„ -ly If 3 dip dip dip 3 dip 

^ij 2 ^^'^ " + (^2 2 dxi dx3 + + 4 dxh ^ 



1 

4^2 

(3.11) 
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Moreover, working out the conservation laws (j2.4ap . the latter are seen to assume 
the form 



1 



1 

2^ 



(3.12) 



Equations ()3.12p allow to clarify the relationship between the conservation laws 
(|2.4|) and the symmetrized Einstein-like equations (|3.11|) . To start with, we recall 
that equations ()3.6p and ()3.7p are equivalent to the relation 



S - -V(ip) + 2V'(ip) = 



(3.13) 



(see [5] for the proof). After that, taking the trace of equations ()3.1ip . we get 



Substituting equation (j3.14p in equation (j3.13p . we obtain 



(3.14) 



R + ^ip.ip'^ - -V.ip^ + ^V{ip) - -V'iip) + -^Sp,SPi = (3.15) 
We rewrite equation (|3.1ip in the form 



^ dip 3 dip dp i^j^ 



(3.16) 



The covariant divergence of (|3.16p yields 



iy'p)Rij + pV'G^j - -RVjp = V*S(,,) + [V'ViVj - VjWij p+ 



95 



95 



(3.17) 
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By definition, the Einstein and the Ricci tensors satisfy V^Gij = and {'V^ip)Rij 
(v^VjVi — ViV^Vj^ ip. Then equation ()3.17p reduces to 



1 ^. 



--iiV,-vP = V*S(,,-) + V* 



9? 



- -T.'^ifi + -7'fh'P 9ij - V{(p)gi 



+ 



V'Vh{S\Uj+S^jU,]+ (3.18) 



On the other hand, making use of equation (|3.15|) it is easily seen that 



1 



-RVjip = V* 



3 3 

+ -iph^^Qij - V{ip)gij 



8p 



ShiS^' (3.19) 



Inserting equations (I3.12p and ()3.19p into ()3.18p . we end up with final equations 



1. 



' ■ 2 



Wh[S\U,+S'',Ui 







(3.20) 



representing the requirement of vanishing for the covariant divergence of the sym- 
metrized Einstein-hke equations, written in the form ()3.16p . As we shall see in the 
next section, equations (j3.20p have to be imposed in order to achieve consistency 
between the assignment of initial data, their preservation in time and the field 
equations. 



4 Bianchi-I cosmological models 

In this section we shall study Bianchi type-I models arising from the above de- 
scribed theory. To this end, we consider a Bianchi type-I metric 

ds^ = df - a\t) dx^ - b\t) dy^ - c\t) dz'^ (4.1) 

The non-trivial Christoffel symbols associated with the metric (|4.ip are 

Setting W = (1,0,0,0), the convective condition (j3.2p implies Soi = 0; from 
this, making use of equations ()4.2p and taking the hypothesis of homogeneity into 
account, it is easily seen that the energy-momentum tensor reduces to 

^ij = (P + P) UiUj -pQij (4.3) 
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as well as that the non-diagonal part of the Einstein- like equations (|3.1ip , evaluated 
for the metric (|4.1|) . yields the equations 



5l2 - - T 



(4.4a) 



Sl3[ 



(4.4b) 



'23 



(4.4c) 



Neglecting the trivial case (zero spin and total isotropy), equations ()4.4p admit 
three kinds of possible solutions: the first one is given by Sij = (zero spin); the 
second one is given by ^ = ^ = ^ (total isotropy); the third one is characterized by 
a partial isotropy and by a spin aligned along a spatial axis. It is clear that the only 
significant solutions are represented by the third case. Indeed, on the one hand in 
the first case we would have an anisotropic universe without spin density, thus the 
anisotropy would be unjustified; on the other hand, in the second and opposite 
case, we would have discrepancy between the geometrical isotropy of the universe 
and the non- vanishing spin density generating a spin vector with non-zero spatial 
components. Therefore, in the subsequent discussion, we shall focus on solutions 
of the third kind only. For instance, supposing the spin aligned along the z-axis, 
a Bianchi type-I metric consistent with the stated assumption is then of the form 
(|4.ip with a{t) = b{t), while 5i2 = —S21 are the only non- vanishing components 
of the spin tensor. In such a circumstance, the diagonal part of the Einstein-like 
equations gives rise to the equations 



+ 2- 



a c 



a c 



-ip -Ifif- 



V{ip) 



1 

V 



(4.5a) 



a c ac 
- + - + — 
a c ac 



a 

2- + 
a \ a 



n 1 
-- + — 



n 1 



a 6 2 I 'T 
a 4 \ T 



if if 
, T := a 

laws (j3.8p and (j3.12p give rise to the equations 



Vif) 



Viif) 



+ ^S^ (4.5b) 



1 



+ _52 (4.5c) 



where we have put r := a^c and S'^ := ^S'^'^Spq [15j. In this case, the conservation 



p+{p + p)- = 

T 



5 + 5- = 
r 



(4.6a) 
(4.6b) 
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Supposing an equation of state of the kind p = Ap (0 < A < 1), the latter admit 
general solutions 

P = Por-^''-'\ S = ^ (4.7) 

r 

po and 5o being suitable integration constants. Moreover, under the stated hy- 
potheses on the metric tensor and the spin tensor, it is an easy matter to verify 
that equations (j3.20p are trivial identities automatically satisfied. 
Concerning the Einstein- like equations (|4.5p , we can subtract equation ()4.5bp from 
equation ()4.5cp obtaining 

'^^i; =0 4.8 

at \a c J \a c J \a c J 

from which we derive 

(4.9) 

c 

D and X being suitable constants, which gives the evolution equation for the 
ratio of the two axes that are not constantly proportional, and therefore it can be 
considered as the evolution equation for the shape of the universe; from equation 
(j4.9p the two explicit expressions follow 

a = DirseVs J -pr ) (4.10a) 

c = Z)~lr3eV 3 J -pr) (4.10b) 

Also, multiplying ()4.5ap by 3 and adding the result to the summation of two times 
equation (I4.5bp with ()4.5cp . we get the final equation for r 

2^ + 3^ = 3^-3^-5^^- (4.11) 

which can be considered as the evolution equation of the volume of the universe. 
Here, it is worth noticing that in equation (j4.1ip (as well as in equation (j4.8p ) 
there is no longer the presence of spin contributions. Another relevant remark is 
that equation ()4.5ap plays the role of a constraint on the initial data: thus for 
consistency we have to check that, if satisfied initially, this constraint is preserved 
in time. To see this point, we first observe that the Einstein-like equations ()3.1ip 
(or (|3.16p ). and thus also (|4.5p . can be written in the equivalent form 

R^i = Tij - Ifg.j (4.12) 



where 



~ _ 1 1 / 3 dip dip \7 d^ dip dip 

:- -^(,,) + -2 ^"2 9^^ + "f^^d^i + Adx^d^^ 

jh dip T/"/',„^„ 1 I ^ ^ ch TT I 1+1 TT \ ^ V7 / ch TT ch 



(4.13) 
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denotes the effective energy-momentum tensor appearing on the right hand side 
of equations while T is its trace. It is then a straightforward matter to 

verify that equations ()4.8p and ()4.1ip can be equivalently obtained by suitably 
combining the space-space equations of the set (|4.12|) : more in detail, equation 
(|4.8p is obtained by subtracting form one another the two distinct space-space 
equations of (j4.12p , while equation (j4.1ip is obtained adding together all the space- 
space equations of ()4.12p . As a consequence, we have that solving equations ()4.8p 
and (|4.1ip amounts to solve all the space-space equations of the set (j4.12p . In 
addition to this, as it is proved in Section 3, the conservation laws ()2.4p together 
with equations (j3.20p automatically imply the vanishing of the four-divergence of 
the Einstein-like equations in the form ()3.16p . The two just mentioned facts allow 
to apply to the present case a result by Bruhat (see [16], Theorem 4.1, pag. 150) 
which ensures that the constraint (|4.5ap is actually satisfied for all time. 



5 From the Jordan to the Einstein frame 

In this section we shall show that, by passing from the Jordan frame {gij) to the 
Einstein frame {gij) through the conformal transformation gij := ipgij (if (/? > 0) or 
gij := —ipgij (if < 0), we can obtain a final differential equation for the volume- 
scale of the universe which is simpler than equation ()4.1ip and, in general, can be 
quadratically integrable. The new equation arising in the Einstein frame allows a 
qualitative analysis of the corresponding solutions even when it is difficult to solve 
analytically. 

To achieve our aim, we have to perform two distinct steps: the first one consists 
in expressing all field equations ()3.8p . ()3.1ip . (|3.12p and ()3.20p in terms of the 
conformal metric gij, its Levi-Civita connection T-^ and the four velocity C7* 
{gijWU^ = I); the second one is evaluating all the so obtained field equations for 
a Bianchi type-I solution of the kind 

= dt^ - a^{t) {dx^ + dy^) - c'^it) dz^ (5.1) 

and for a spin tensor whose unique non- vanishing components are again S12 = —S21 
(spin aligned along the z-axis). From now on, for simplicity, we suppose that ip > 0. 
Introducing the conformal metric gij = pgij, it is an easy matter to verify that 



Gij — Gij ^2 



3-3 

+ f^jPi + -^PhP^gij - <p{^^^h)gij 



(5.2) 



where Gij and Gij are the Einstein tensors associated with the metrics gij and gij 
respectively. We can then write the symmetrized Einstein-like equations (|3.1ip in 
the simpler form 



Gij = - ^V{p)gij + ^ (nS'^jU, + fhS\Uj ) + 

1 / N 1 1 

— ' TT ah ' " ^ - 



(5.3) 
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Equations (j5.3p are the key point of our discussion. To start with, making use of 
the relation 

f > = f > + ^^,6^ - ^^,9''9., + ^^.S^ (5.4) 

Unking the Levi-Civita connections associated with the metric tensors gij and gij 
respectively, and taking the identities W = sJ^lJ^ {Ui = -^Ui) into account, it is 
easily seen that under the stated hypotheses the energy-momentum tensor ()3.10p 
assumes the form 

1 - - p _ 
^ij ■■= -{p + p)UiUj Qij (5.5) 

if LP 

Inserting the content of equations ()5.ip . ()5.4p and ()5.5p into equations (|5.3p . we 
get 

2 

-\ +2^^_ = A^-\V{^)-^S' (5.6a) 



and 

a c a c {p'^ [f-^ 4(/9'^ 



2r+(Z) =_J__^v{^) + -^S' (5.6c) 



In a similar way, from the conservation laws (jS.Sp and (j3.12p we obtain the equa- 
tions 

iS r 2 (/p 

^ + (1 + A);-|^(1 + A) = (5.8) 

p T lip 



and 



Setting := ib(/3 2^ one has —\^ = % and then we can integrate equation (|5.7p 
and (|5.8p as 

5=^, p = Mre)-(i+^) (5.9) 



where So and /Oq are integration constants. Concerning equations (|5.6p . we can 
proceed as in the Jordan frame (see section 4), so obtaining the relations 

^ = i^exp [x \\dt\ (5.10) 



and 

2^ = 34-34--^y(v.) (5.11) 

r ip^ ip^ ipi^ 

where D and X denote suitable integration constants. We stress that, as in equa- 
tion (|4.1ip . also in equation (|5.1ip there are no spin contributions. At this point, we 
notice that if we could express the quantities p, and ip (or, equivalently, = (/7~2 ) 
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as functions of r and insert the so obtained results into equation (jS.lip . we would 
get a final equation for r of the kind 



fir) 



(5.12) 



which is quadratically integrable. The second equation of (|5.9p together with 

2 

Q~3 = = f'[R(^p)) constitute a set of two relations involving the three variables 
f,peQ. In principle, it is then possible to express the last two variables as suitable 
functions of r alone, thus obtaining from (|5.1ip a final quadratically integrable 
equation of the kind (|5.12p . 

As a last remark, we note that equation ()5.6ap has the nature of a constraint 
on the initial data. As made in the Jordan frame, we have again to check that 
such a constraint is satisfied for all time after the initial instant. In connection 
with this, we now prove that the conservation laws (|2.4p together with equations 
()3.20p ensure that the four-divergence (with respect to the Levi-Civita connection 
associated with the conformal metric gij ) of the effective energy-momentum tensor 



I. 



-s\u, 



S\Ui 



2ip 



■^-^ S^*^ SpgUiUj + -^^S^'^Spggij 



(5.13) 



vanishes. To see this point, indicating by Vj the Levi-Civita covariant derivative 
associated with the conformal metric gij, we have 



•P 



dip 

2ip yOx^"'' ' 'dx^^^ 



Q^u9 9sj 



(5.14) 



We have separately 



if-' 



1 



1 



-^^S^'^SpqUiUj + -^^S^'^Spqgij 



J 



^^S^'^ SpqUiUj + ^^S^'^ Spqgij 



+ 



2(/pa 



if-' 



-S%Uj - S^:Ui + J—SP'^SpqUiUj - ^SPiSpqg.j 



(5.15a) 



12 



2^2 



if 



'Pi 



4^ 



V{p) 



<Pi 



hk 



■ShkS' 
(5.15b) 



1 



1 



+nS%Uq) + —V,, (-S\Uq - S\Ui 



2ip 



1 



■S^^ShkUiUq + -^-^S^^Shkdiq 



^S^'^SpqUiUj + ^^^S^'^SpqQij 



(5.15c) 



Collecting equations ()5.15ap . (|5.15bp and ()5.15cp . we end up with 



z 99 



^ \ -^^ShkS^^ 1 + 



(5.16) 



4ip 



Inserting the content of equations (j3.12p and (jS.lSp into (j5.16p . we conclude that 
(holding the conservation laws ()2.4p ) the requirement V-^T,j = is mathematically 
equivalent to the conditions 



1. 



1. 







(5.17) 



identical to the already obtained equations ()3.20p . The last step (left to the reader) 
is to verify that, for the conformal metric (|5.2p and the given spin tensor, equations 
p.20p are again trivial identities. 

Repeating the arguments developed in section 4 for the Jordan frame, we can apply 
again the Bruhat's result [16] and thus verify that the Hamiltonian constraint 
(|5.6ap is actually preserved in time. 
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6 Examples 



2 



Example 1: f{R) = model. The first example we consider is given by the 

2 

model f{R) = Rs , discussed in the Einstein frame. In this case, the associated 
potential ()3.7p is expressed as 

V{ip) = -^^-^ (6.1) 
while the scalar field ()3.6p is given by 



1 



^=^P- (6.2) 

We assume an equation of state for the fluid of the kind p = ^p. From equation 
()5.9p we have then 

p = po{f Q)-l (6.3) 

3 

Equations (j6.2p and (j6.3p together with the identity Q = (f~2 yield the relation 



fe = ( ) po^ri (6.4) 



1 4 



from which we get the final expressions for the scalar field 



4 

9 2 10 



V>= { Po ^r27 (6.5) 



and the energy density 



10 
9 



4 20 



Inserting the content of equations (|6.ip . (|6.5p and (|6.6p into equation (jS.lip . we 
end up with the final dynamical equation 
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f = Bf-27 (6.7) 

where B is a positive constant (depending on pQ.) Integrating a first time, we get 
the relation 

1 27 

-f^ = -Bf^+E (6.8) 



E being an integration constant. For values E < 0, equations (j6.7p and (j6.8p 
imply the existence of a minimum for the volume-scale f, so giving rise to a 
universe without singularity. Moreover, by plotting the functions r and j, it is 
seen that the volume-scale f is infinite of order 1 < n < 2; due to equation ()5.10p 
this fact implies that the universe isotropizes at large t. Also, at large t, we have 
^ ^ f 27 and then the universe expands also with respect to the Jordan 

frame gij = ^Qij- In connection with this, we notice that the integral J^* '^^ 
diverges for t — )• +oo, thus the time reparametrization di = -^dt (making the 
Jordan frame a Bianchi-I metric) is well defined. 
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Example 2: f{R) = model. We discuss a second example in order to show 
how solving the field equations in the Jordan or in the Einstein frame are in general 
two non-equivalent procedures yielding different results. To this end, we consider 
the model f{R) = R'^ coupled with a dust fluid (A = 0). In such a case, the 
potential ()3.7p is expressed as 

Viv) = ' (6.9) 

In the Jordan frame, using the first equation of (j4.7|) . it is easily seen that the 
scalar field ()3.6p is given by 

(^ = 2tp|r"3 (6.10) 

Inserting equations (j6.9p and (j6.10p into equation (j4.1ip . we end up with the final 
equation 

Jf = +Ai! + c4rf (6.11) 

9 

where C = 7(2)~4 is a positive constant. Equation (j6.1ip admits solutions of the 
form 

T{t) = {At + Bf (6.12) 

with A and B constants. Indeed, substituting equation (j6.12p into equation (j6.1ip . 
we have 

(^2A^ -Cp^^ {At + Bf = (6.13) 

which is always satisfied by suitably choosing the constant po as function of the con- 
stant A or vice versa. From equations (|6.1Up and (|6.12p we have np ~ {At + B)"^; 
due to equation ()4.9p and ()4.10p . at large t solutions ()6.12p are seen to isotropize 

8 

as well as they undergo an accelerated expansion, being a{t) w c{t) w {At + B)^. 
It has not escaped our attention that a model of dust may not necessarily be the 
most adequate to describe the universe in its early stage; for such a situation, 
where quantum effects might arise, the coupling to the Dirac field seems to be 
more appropriate. In this respect, it is worth noticing that in the case of the 
Dirac field ^p it is possible to find a solution analogous to (|6.12p : indeed in such a 
circumstance the Dirac field density '^l>^|J is seen to play the same role of the energy 
density p, and the eventual dynamical equation has the same form of (|6.1ip . We 
also stress that in the presence of the Dirac field, the existence of the spin density 
makes the choice of the initial anisotropy unavoidable; inasmuch as the universe 
has to evolve according to an initial infiationary expansion together with a phase of 
isotropization, the solution ()6.12p appears to give a plausible model. For a better 
understanding of this point, we refer the reader to reference [11| . 
Now we can study the same model in the Einstein frame. This time we have 

= 2tp| (fe)"3 (6.14) 
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3 

Making use of the definition = ip'^ 2 , from equation ()6.14p we derive the relation 

e = 2i5(^)' (6.15) 

from which we get 

34 = 3(2)1^0^"^ (6-16a) 

- \V{^) = -9(2)V^t-^ (6.16b) 

For values of po sufficiently small, inserting the content of equations (|6.16p into 
equation (|5.11|) we obtain a final equation of the form 

f = (6.17) 

C being a positive constant (dependent on po). Equation ()6.17p can be easily 
integrated as 

A{t + By + ^ (6.18) 

where the constants ^ > and E have to be determined consistently with equation 
(j5.6ap . From equation (j6.18p we have > ^: we see that in the Einstein frame 
the model f{R) = R'^ coupled to dust can give rise to a universe with no initial 
singularity but undergoing a slow expansion without isotropization. 

3 

Example 3: f{R) = R^ model. As third and last example we discuss the 
model f{R) = R'^ coupled with a fluid having equation of state of the kind p = f . 
To start with, it is easily seen that the potential ()3.7p assumes here the form 

Viv) = (6.19) 

Moreover, from the trace f|2.2p of the Einstein-Hke equations we deduce the relation 

2 

R = and then from equation f|3.6p we have 

if = -ph (6.20) 

Inserting the content of equations (j6.19p and (j6.20p into the dynamical equations 

(|4.1ip and (jS.lip . a direct calculation shows that in both cases the contribution 

due to density and pressure elides that due to the potential V{ip). The conclusion 

i 1 

follows that in the Jordan frame (taking the relation (/? = |poT~2 into account) 
we end up with the dynamical equation 

1 ^2 

(6.21) 



2 r 
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while in the Einstein frame we have 

f = (6.22) 

The general (positive) solution of (j6.21|) is given by r(t) = (At + B), while the 
solution of ()6.22p is of course f{t) = At + B (A and B being integration constants). 
We then see that both solutions in the Jordan and in the Einstein frame have 
singularity and do not isotropize. 

7 Conclusion 

In the present paper we have studied anisotropic cosmological models of the kind 
Bianchi type-I arising within the framework of /(i?)-gravity with torsion coupled 
with Weyssenhoff spin fluids. Following the lines traced in [llj for the coupling 
with Dirac spinors, we have derived and discussed the field equations in both the 
Jordan and the Einstein frame, showing that in the latter the Einstein-like equa- 
tions can be put in a quadratically integrable form. This fact can turn out to be 
useful for a qualitative analysis of the solutions or numerical integration even in 
the case analytical integration results hard. 

As it has been noticed in previous works [3 [H [3 [8] , in /(i?)-gravity with torsion 
there are two possible sources for torsion: the spin density and the non-linearity of 
the gravitational Lagrangian function f{R). In general, the presence of non-zero 
spin forces us to take into account anisotropic models, since in such a case torsion 
does not satisfy the cosmological principle [T7]; on the contrary, when spin is zero, 
the condition of homogeneity gives rise to a torsion consistent with the hypothesis 
of spatial isotropy, in such a circumstance the torsion being generated by a time- 
vector with zero spatial components. 

Anyway, dealing with Weyssenhoff spin fluids and Bianchi type-I space-times, we 
have shown that the only significant model is that characterized by spin aligned 
along a spatial axis and partial isotropy (as it happens in the already known case 
f{R) = R). Moreover, we have proved that for such models the role played by 
the spin is circumscribed to the assignment of the initial data; indeed, in both 
the Jordan and Einstein frames, spin contributions are seen to disappear from 
the dynamical equations (j4.8p and (j4.1ip or (|5.1ip . while spin is present in the 
Hamiltonian constraints ()4.5ap and ()5.6ap : the only torsional contribution to the 
dynamics seems then given by the terms arising from the non-linearity of the func- 
tion f{R); such terms are present in an identical form also in the case of zero spin. 
In addition to this, we have studied a few examples, by picking functions of the 
type f{R) = R"' in different cases, and we have seen that the resulting theory 
seems to be quite sensitive to the specific choice of the model. 
Finally, we have derived general conditions assuring the vanishing of the quadridi- 
vergence of the Einstein- like equations in both the frames; such conditions turn out 
to be identically satisfied by the models here considered and allow the preservation 
in time of the Hamiltonian constraint (momentum constraints are trivial identi- 
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ties); therefore we have proved the correct formulation of the associated Cauchy 
problem in both the Jordan and Einstein frames. 
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